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基础 研究 
基于 流 形 学 习 的 血管 内 超声 图 像 序列 关键 帧 的 提取 太 应 用 
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摘要 :目的 以 流 形 学 习 为 基础 ,提出 一 种 基于 血管 内 超声 GVUS) 图 像 序列 的 关键 帧 门 控 方法 ,抑制 IVUS 序列 纵 切 方向 上 的 运 
动 伪 影 。 方 法 应 用 流 形 学 习 方 法 中 的 拉 普 拉 斯 特征 映射 算法 ,将 高 维 IVUS 图 像 序列 降 到 低 维 流 形 中 ,利用 低 维 特征 向 量 , 构 
建 一 个 距离 函数 来 反映 心脏 运动 规律 ,将 IVUS 图 像 分 为 心脏 舒张 未 期 和 非 心 脏 舒张 未 期 两 类 ,从 而 提取 关键 帧 ,组 成 门 控 序 
列 。 结 果 临床 采集 13 组 IVUS 序 列 (图 像 915+142 帧 ,血管 长 度 15.24+2.37 mm) ,计算 门 控 前 后 图 像 序列 的 血管 容积 、 管 腔 容 
积 和 平均 斑 块 负荷 。 统 计 实验 结果 ,表明 门 控 序列 血管 容积 、 管 腔 容 积 显 著 小 于 原始 序列 , 门 控 前 后 序列 的 平均 斑 块 负荷 差异 
性 不 显著 ,满足 临床 诊断 要 求 。 血 管 面积 方差 和 管 腔 面 积 的 方差 显著 小 于 原始 序列 ,表明 门 控 序列 较 原始 序列 稳定 。 在 IVUS 
图 像 序 列 的 纵 切 图 像 上 , 门 控 序列 减少 了 锯齿 形状 的 运动 伪 影 ,与 原始 序列 形状 一 致 , 且 具 有 良好 的 连续 性 。 并 将 本 文 方法 与 
已 有 的 提取 门 控 序列 方法 进行 对 比 。 结 论 本 文 方法 算法 简单 稳定 ,抑制 了 IVUS 图像 序列 的 纵向 运动 伪 影 。 

关键 字 : 血 管内 超声 成 像 ; 门 控 ; 关 键 帧 ; 流 形 学 习 


Key frames extraction and application in intravascular ultrasound pullback sequences 


based on manifold learning 


MAO Haiqun', YANG Feng’, LIN Mudan’, HUANG Zheng’, CUI Kai, WANG Xinxin' 
‘College of Biomedical Engineering, Southern Medical University, "Department of Cardiology, Southern Medical University, Nanfang 
Hospital, Guangzhou 510515, China 


Abstract: Objective We propose an image-based key frames gating method for intravascular ultrasound (IVUS) sequence 
based on manifold learning to reduce motion artifacts in IVUS longitudinal cuts. Methods We achieved the gating with 
Laplacian eigenmaps, a manifold learning technique, to determine the low-dimensional manifold embedded in the 
high-dimensional image space. A distance function was constructed by the low-dimensional feature vectors to reflect the heart 
movement. The IVUS images were classified as end-diastolic and non-end-diastolic based on the distance function, and the 
IVUS images collected in end-diastolic stage constitutes the key frames gating sequences. Result We tested the algorithm on 13 
in vivo clinical IVUS sequences (images 915-142 frames, coronary segments length 15.24 2.37 mm) to calculate the vessel 
volume, lumen volume, and the mean plaque burden of the original and gated sequences. Statistical results showed that both 
the vessel volume and lumen volume measured from the gated sequences were significantly smaller than the original ones, 
indicating that the gated sequences were more stable; the mean plaque burden was comparable between the original and gated 
sequences to meet the need in clinical diagnosis and treatment. In the longitudinal views, the gated sequences had less saw 
tooth shape than the original ones with a similar trend and a good continuity. We also compared our method with an existing 
gating method. Conclusion The proposed algorithm is simple and robust, and the gating sequences can effectively reduce 
motion artifacts in IVUS longitudinal cuts. 

Key words: intravascular ultrasound imaging; gating; key frames; manifold learning 


TEE S KP LC Co ML AER TIS] ESI EROR 影 (CAG) 非 常 重要 的 辅助 工具 。 然 而 ,由 心脏 运动 所 引 
血管 内 超声 (intravascular ultrasound, TVUS) 是 一 种 基 起 的 运动 伪 影 给 IVUS 序列 血管 面积 .血管 容积 . 管 腔 面 
于 导管 的 成 像 技术 ,实时 提供 冠 脉 血管 的 截面 图 像 "。 H 管 腔 容积 管 腔 最 大 最 小 直径 \ 管 腔 面积 狭 罕 率 、 斑 块 
图 像 分 辨 率 高 ,能 准确 获取 冠状 动脉 的 管 壁 厚度 . 管 腔 负荷 等 诊断 指标 的 定量 分 析 华 来 很 大 的 困难 ,影响 临床 
形态 及 动脉 粥 样 硬化 斑 块 成 分 等 信息 。 在 辅助 诊断 冠 “”” 冠 心病 的 诊断 和 治疗 "> 。 


脉 粥 样 硬化 方面 ,IVUS 已 经 发 展 成 为 传统 冠状 动脉 造 运动 伪 影 在 IVUS 图 像 序列 中 表现 非常 显著 ,在 纵 
切 图 上 ,IVUS 图 像 序 列 呈 锯齿 状 分 布 (如 图 6 上 部 所 
收 稿 日 期 :2014-10-05 I)o ECG 门 控 法 是 传统 的 抑制 运动 伪 影 的 方法 ,分 为 
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KR ECG 信号 和 IVUS 图 像 ,医生 对 照 ECG 信号 从 
IVUS 序 列 中 标记 出 门 控 帧 ,介入 时 间 减 少 ,缺点 是 难以 
选择 最 佳 采 样 点 ,主观 性 强 。 临 床上 难以 同时 获得 
ECG 信 和 号 和 IVUS 信 号 “。 因 此 ,近年 来 ,人 研究 者 们 致力 
于 利用 IVUS 图 像 序列 本 号 来 寻找 门 控 序 列 , 即 基于 图 
像 的 门 控 法 “2 。 基 于 图 像 的 门 控 法 完全 利用 图 像 处 理 
技术 ,提取 反映 TVUS 图 像 的 显著 性 特征 量 作 为 参考 点 ， 
参考 点 随 着 图 像 序列 变化 ,形成 一 个 类 似 ECG 的 时 间 
变化 信号 。 据 此 将 IVUS 图 像 分 为 心脏 舒张 未 期 图 像 
和 非 心脏 舒张 未 期 图 像 两 类 ,由 关键 帧 组 成 门 探 序列 。 
所 谓 的 关键 帧 定义 为 心脏 舒张 末期 采集 到 的 IVUS 图 
像 。 在 心动 周期 中 ,舒张 未 期 心室 进入 等 容 收缩 期 , 心 
脏 运动 较 慢 ,超声 探头 和 管 腔 的 相对 位 移 小 ,此 时 的 
IVUS 岁 像 具有 很 好 的 稳定 性 和 一 致 性 "2 。 目 前 已 有 
的 基于 图 像 的 门 控 方 法 有 :以 像素 点 局 部 灰 度 均值 " , 感 
兴趣 区 域 (ROI, 人工 设 定 ) 的 平均 灰 度 和 绝对 灰 度 差 "， 
管 腔 轮 廓 线 ” 等 局 部 信息 作为 特征 的 方法 。 文 献 [6] 提 
出 的 方法 对 带 通 滤波 带 参 数 敏感 ,难以 提取 心律 不 齐 患 
者 的 IVUS 序 列 门 控 帧 ,文献 [7-8] 提 出 的 方法 都 需要 精 
准 的 边缘 检测 方法 ,然而 目前 尚未 存在 有 效 的 IVUS 图 
像 分 割 方法 ,降低 了 特征 精度 。 文 献 [9] 提 出 一 种 基于 
图 像 的 “理解 (Intelligate)” 的 门 控 方法 ,选择 图 像 互 相关 
系数 , 空 频 域 特征 等 全 局 特征 作为 特征 信息 ,缺点 是 难 
以 准确 选择 主要 特征 。 为 了 降低 对 特征 提取 与 选择 的 
要 求 ,研究 者 利用 同 序列 中 不 同 帧 之 间 的 非 相似 性 程度 
作为 特征 信息 ”” ,构造 差异 矩阵 ,采用 动态 规划 
(dynamic programming, DP) ,在 差异 矩阵 中 搜索 一 条 
累计 差异 值 最 小 的 路 径 ,提取 出 心脏 舒张 未 期 门 控 帧 ， 
算法 复杂 ,计算 量 大 。 

流 形 学 习 是 一 种 非 线性 降 维 方法 ,将 高 维 非 线性 数 
据 映 射 到 低 维 流 形 , 同 时 保持 了 数据 集 的 固有 属性 ,在 
图 像 处 理 中 得 到 广泛 应 用 "在 门 控 方面 , 流 形 学 习 已 
用 于 IVUS 心 脏 门 控 外 .4D 超声 “ ”MRI 图 像 的 呼吸 门 
探 ”。 本 文 提出 一 种 基于 流 形 学 习 的 关键 帧 门 控 方 法 ， 
在 流 形 学 习 的 指导 下 ,将 高 维 非 线 性 的 IVUS 数 据 集 降 
到 低 维 流 形 中 ,利用 低 维 流 形 中 的 特征 向 量 ,构造 一 个 
特征 明显 的 距离 函数 ,来 反映 心脏 运动 规律 ,检索 距离 
函数 局 部 最 小 值 ,将 关键 帧 从 IVUS 图 像 中 提取 出 来 ,组 
成 门 控 序列 。 通 过 对 门 控 前 后 序列 的 血管 面积 ,血管 容 
积 , 管 腔 面 积 . 管 腔 容积 、 途 块 负 答 等 量化 指标 的 定量 比 
较 分 析 , 以 及 门 控 前 后 序列 纵 切 图 分 布 和 与 已 有 的 方法 
比较 来 评价 本 文 算法 。 


1 方法 
TEIVUS 图 像 序列 采集 过 程 中 ,心脏 周期 性 搏动 引 
起 超声 探 尖 与 血管 之 间 的 相对 位 移 、 短 轴 方 向 上 的 平面 


振荡 、 血 管 尺寸 连续 变化 、 冠 脉 血管 部 分 位 置 重 采样 或 
未 采样 等 ,在 IVUS 纵 切 图 像 上 产生 锯齿 形状 的 运动 伪 
影 。 但 这 些 变化 都 是 缓慢 平滑 的 , 且 同 一 相位 处 的 
IVUS 图 像 具有 很 高 的 相似 性 。 因 此 ,IVUS 图 像 序 列 处 
于 一 个 环绕 空间 (ambient space) 的 低 维 平滑 流 形 上 ， 
一 帧 图 像 代 表 高 维 空间 中 的 一 个 样本 点 。IVUS 图像 序 
列 具 有 数据 量 大 、 维 数 高 . 非 线 性 的 特点 ,一 幅 图 像 中 包 
含 成 千 上 万 的 图 像 特征 。 心 脏 运动 只 需 用 儿 个 基本 参 
数 描述 就 能 实现 ,有 效 选 择 反 映 心脏 运动 的 特征 参数 是 
本 文 的 重要 内 容 。 本 文采 用 拉 普 拉 斯 特征 映射 算法 
(laplacian eigenmap,LE) 对 IVUS 序列 图 像 进行 降 维 ， 
LE 算法 保持 了 数据 流 形 内 在 的 几何 特性 ,求解 过 程 可 
转化 为 对 稀 跑 矩阵 的 特征 分 解 问题 ,计算 简单 。 
1.1 拉 普 拉 斯 特征 映射 算法 

拉 普 拉 斯 特征 映射 算法 的 基本 思想 是 利用 图 的 谱 
理论 ,寻找 某 种 意义 上 最 优 地 保持 局 部 邻 域 信息 的 低 维 
表达 ,使 得 在 高 维 空间 中 距离 近 的 点 在 映射 到 低 维 流 形 
后 它们 之 间 的 距离 仍旧 很 近 "”。 给 定 高 维 空间 R^ 中 的 
一 组 样本 点 X={x,…,xj), 寻找 低 维 空间 RN >> 1) PR 
一 组 点 7={y1,.…,y,}, 用 低 维 y 描 述 高 维 x;。LE 算 法 主要 
由 三 个 步 又 组 成 : 

(1) 选 择 邻 域 ,构建 局 部 邻 域 图 G。 计 算 每 个 样本 
点 x 和 其 余 样本 点 之 间 欧 氏 距 离 ,如 果 样 本 点 x 和 x 是 近 
邻 点 , 则 图 G 有 边 x;x;。 可 采用 以 下 两 种 方法 确定 邻 域 : 

Ok - nr 算法 :如 果 x 是 x 的 & 近 邻 点 或 者 x 是 x 的 k 近 
邻 点 , 则 x 和 x 有 边 连 接 ; 

e - bal 算法 :车 x 与 x 的 距离 小 于 e, 则 点 x 和 xj 有 
边 连 接 , 与 x 相连 接 的 点 都 是 其 近邻 点 。 

(2) 选 择 权 值 W。 在 邻 域 关系 图 G 中 ,设置 权 值 ,有 
以 下 两 各 方法: is 


OD 高 斯 核 方 法 :W=1e | 


,与 近邻 ,其 中 4 


| 0 ， 其 它 
为 一 个 调节 参数 ; ts 
@ 简 单 核 方法 :W=1、 0 T 
0 其 它 


(3) 特 征 映射 ,最 小 均 方 误差 ,计算 样本 的 低 维 具 
入 。 为 保持 数据 的 局 部 近邻 关系 ,最 小 化 误差 国 数 : 


min S W,|y;-y,||. Jh jo 1,2, n, nRa 
ij 


(1) 
在 YDY’ =7 约 束 下 ,公式 (1) 可 转化 为 广义 特征 值 
的 求解 问题 
LY-ADY (2) 


公式 (2) 中 ,L=D-W 为 拉 普 拉 斯 矩阵 ， 
Di = 2 Wo Y ALD 前 :个 不 为 零 的 最 小 广义 特征 值 对 


应 的 广义 特征 向 量 组 成 的 矩阵 。 
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1.2 基于 LE 的 TVUS 图 像 序列 关键 帧 提取 算法 

在 LE 的 指导 下 ,提取 原始 序列 关键 帧 ,主要 过 程 
为 :(1) 坐 标 转换 。 将 直角 坐标 下 的 IVUS 图 像 转换 到 极 
坐标 下 ,去 除非 成 像 区 域 信息 ,减少 计算 量 ;(2)LE 降 
维 。 采 用 LE 算法 ,将 高 维 的 IVUS 数据 集 降 到 低 维 流 
形 中 ,提取 低 维 特征 向 量 ;(3) 距离 函数 构造 。 利 用 选择 
的 特征 向 量 构造 一 个 距离 函数 ,反映 心脏 的 运动 规律 ; 
(4) 提 取 关 键 帧 。 检 索 距离 函数 的 局 部 最 小 值 , 将 关键 
帧 从 IVUS 图 像 中 提取 出 来 ,组 成 门 控 序列 。 
1.2.1 IVUS 坐标 转换 流 形 学 习 的 第 一 步 是 构造 邻 域 
距离 矩阵 ,要 求 对 所 有 的 数据 进行 存储 。 对 于 一 组 像素 
为 512x512, 长 度 为 1024 的 IVUS 图 像 序列 ,所 占用 的 
内 存 空间 是 2 GB, 这 对 计算 机 存储 空间 和 程序 运行 时 
间 的 损耗 都 是 很 大 的 。 首 先 对 图 像 进行 坐标 转换 ,将 
IVUS 图 像 转 换 到 极 坐 标 下 存储 运算 。 同 一 组 序列 , 极 
坐标 下 占用 的 内 存 为 0.7031 G, 只 占 直角 坐标 系 下 容量 
的 35%。 统 计 直 角 坐 标 系 和 极 坐 标 下 LE 降 维 后 ,其 特 
征 向 量 之 间 的 绝对 误差 仅 是 0.9%。 因 此 ,对 IVUS 网 像 
序列 进行 极 坐 标 转换 后 , 既 能 保持 原 有 的 特质 ,又 能 减 
少 计算 机 负担 。 图 1 为 IVUS 图 像 坐标 转换 图 ,A 为 直 
角 坐 标 系 下 的 IVUS 图 像 ,B 为 对 应 的 极 坐标 图 像 。 


图 1 IVUS 图 像 坐 标 转换 


Fig.1 Coordinate conversion for an IVUS image. 


1.2.2 LE 降 维 在 一 组 包含 多 个 心动 周期 ,长 度 为 n 的 
IVUS 图 像 序列 w={w,.…,w) 中 ,对 每 一 帧 图 像 uj, 希望 寻 
找 一 个 特征 问 量 G 来 描述 反映 uw 的 心脏 运动 规律 信息 。 
利用 流行 学 习 理 论 , 令 原始 高 维 空间 x; =w, 低 维 流 形 空 
间 y;=$;。 高 维 空间 的 维 数 N 等 于 整 幅 图 像 像 素 , 低 维 
流 形 的 维 数 为 :。 这 样 , 流 形 学 习 表 示 为 : 

m —> R' 

u; $i 

在 心脏 周期 中 ,每 一 帧 图 像 用 一 个 列 向 量 表示 , 代 
表 高 维 空间 中 的 一 个 点 , 随 着 心脏 的 运动 ,序列 图 像 沿 
着 高 维 空间 或 低 维 流 形 周期 性 地 变化 。 应 用 流 形 学 习 
降 维 最 基本 的 问题 是 维持 局 部 信息 的 一 致 性 ,使 得 相似 
的 图 像 在 低 维 流 形 中 有 相 邻 的 分 布 。 心 脏 运动 信息 包 
含 在 拉 普 拉 斯 矩阵 志 的 前 :个 非 零 的 最 小 广义 特征 值 所 
对 应 的 特征 向 量 中 。 


4i-l,..en (3) 


邻 域 选择 是 LE 算法 的 关键 内 容 ,目前 已 有 的 多 数 
邻 域 选择 算法 只 适应 于 数据 点 分 布 均匀 的 流 形 ,对 于 结 
构 复 杂 的 流 形 不 能 自 适应 的 选择 合理 的 邻 域 ,主要 依据 
经 验 来 选择 邻 域 大 小 。 本 文采 用 k- nn 分 域 算 法 ,选择 令 P 
TK > neo,noa 为 一 个 心动 周期 内 的 IVUS 图 像 帧 数 。 
IVUS 图 像 序列 的 欧 氏 距离 矩阵 呈 一 定 的 周期 性 ,这 个 
周期 就 是 心脏 运动 的 周期 ni。 图 2 为 同一 组 IVUS 图 
像 序列 的 前 200 帧 (A) ,前 1024 帧 (B) 图 像 的 欧 氏 距离 
和 矩阵 , 灰 度 值 越 低 ,表明 图 像 的 距离 越 近 。 由 患者 的 心 
率 f 来 确定 nw,f 选 择 过 程 如 下 :(1) PS IK EGER E 
函数 为 EF:( 图 2 中 第 ; 行 或 者 第 ; 列 ) ,对 E 做 传 里 叶 变 换 ， 
的 幅度 谱 在 正常 心率 范围 内 (50-140 beats/min) 的 峰 
值 所 对 应 的 频率 就 是 图 像 w 估 计 得 到 的 心率 方 (2) 所 有 
图 像 估 计 得 到 的 心率 均值 即 为 心率 估计 值 : 


f= LX fo 这样, 一 个 心动 周期 内 的 IVUS 图 像 帧 数 


为 :Noyce = round ( fsamlfe) , py 每 分 钟 采集 IVUS 图 像 
的 帧 数 。 


Oo ë o 
图 2 同一 组 IVUS 序 列 前 200 帧 (4) , 1024 toi (B) EE EB EE RR 
离 矩 阵 


Fig.2 Visualization of the euclidean distance matrix from the 
first 200 frames (A) and 1024 frames (B) of an IVUS set. 


LE 算法 的 另 一 个 主要 内 容 是 低 维 流 形 的 维 数 ;的 
确定 ,目前 还 没有 固定 的 方法 指导 低 维 维 数 的 选择 。 本 
文 实验 表明 ,对 IVUS 图 像 序列 LE 降 维 , 当 t 大 于 4 时 ,对 
中 离 函 数 影响 不 大 ,因此 本 文选 择 内 药 维 数 1= 4。 

1.2.3 构造 距离 函数 利用 LE 算法 对 IVUS 图 像 序列 降 
维 ,IVUS 序列 在 低 维 流 形 中 的 分 布 主要 有 两 种 情况 ,如 
图 3 所 示 。 图 3 为 两 组 IVUS 序列 降低 到 2D .3D 空 间 的 
分 布 图 绿色 圆 实 点 "是 图 像 在 低 维 流 形 中 的 分 布 ， 黑 
色 六 角 星 符 ” 是 文献 [6] 方 法 寻找 到 的 关键 帧 的 分 布 。 
图 3 上 部 虚 框 内 对 应 的 IVUS 图 像 序列 血管 结构 变化 
大 ,序列 进展 方向 上 出 现 斑 块 且 血管 面积 变化 大 , 低 维 
流 形 中 关键 帧 分 布 离散 ;图 3 下 部 黑色 虚 框 内 对 应 的 
IVUS 图像 序列 血管 结构 变化 不 明显 ,无 广 块 ,序列 进展 
方向 上 血管 面积 变化 小 , 低 维 流 形 中 关键 帧 分 布 集中 。 
统计 13 组 临床 TVUS 图 像 序列 ,其 中 10 组 序列 低 维 分 
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布 属 第 一 种 情况 , 另 3 组 序列 低 维 分 布 属 第 二 种 情况 。 
可 见 , 在 低 维 流 形 中 , 仅 用 样本 点 分 布 的 聚 类 性 质 很 难 
寻找 关键 帧 。 文 献 115] 将 低 维 流 形 中 前 后 帧 欧式 距离 
函数 的 局 部 极 大 值 作为 关键 帧 所 处 位 置 , 但 欧式 距离 函 


XS pus (2) 向 前 搜索 局 部 最 小 值 。 以 py 为 搜索 起 
点 ,nwa 为 搜索 步 长 。round(0.4xnoa) 为 窗 宽 , 问 前 搜 
索 上 一 个 局 部 最 小 值 点 pwuwr-ie PEAD siari- 为 新 的 起 点 ， 
寻找 py 的 上 一 个 局 部 最 小 值 点 ps,,.s, 重 复 此 过 程 ， 


数 特征 不 明显 ,容易 漏 检 或 误 检 关键 巾 。 根 据 心脏 运动 
的 特点 ,心脏 舒张 未 期 的 [VUS 图像 稳定 , 相 邻 帧 之 间 的 
差异 性 小 。 本 文 构造 距离 函数 来 提取 关键 帧 ,距离 函数 
定义 为 : 


直到 向 前 搜索 完毕 ;(3) 向 后 搜索 局 部 最 小 值 。 类 似 (2) 
的 方法 ,以 pw 为 起 点 ,向 后 搜索 局 部 最 小 值 点 ;(4) 所 
有 局 部 最 小 值 点 位 置 上 的 图 像 构 成 关键 帧 ,形成 门 控 
序列 。 


dÙ) = wY Md sisi Lat (4) 图 4 中 ,“ 星 号 ”点 为 检索 到 的 距离 函数 局 部 最 小 
e 值 , 即 关键 帧 图 像 所 在 位 置 。 
0.08 = T? 
i p 2 实验 结果 的 讨论 和 分 析 
sd : á M 临床 实验 数据 来 源 于 南方 
£5 0.04 医院 7 名 患者 (4 男 3 女 ,年龄 
0.02 Ab d PC 0.02 e. p 
i y". ^ ma 0 T om 62+9 岁 ,图 像 915+142 帧 ,血管 
RR ge T e K HE 2524 +2.37 mm, 像素 
0 iy Ki "T4: . e Pu 
"T f v i sid. oam. 006 23 p^ 512x512,256 lr Jk FE ) , 13 £H 
| H | A ON £p 05 IVUS 图 像 序列 。 数 据 采 集 设 
-0. " » * » è " t ` y S a = 
-0.06 -0.04 -0.02 0 0.02 0.04 0.06 0.08 pc s B 备 是 是 波士顿 科学 公司 的 
iLab™ Ultrasound Imaging 
0.08 0.06 -— NS ENAT 
xin 004 System, 超 声 探头 为 机 械 式 旋转 


A 
0.04 "cs s 0 
eden c 
j 0.02 
M 0.04 
-005 -0.03 -0.01 0 0.01 0.03 0.04 04 0106-00 


3 两 组 IVUS 序 列 2D、3D 低 维 流 形 分 布 
Fig.3 2D and 3D manifold of two IVUS sequences. 


公式 (4) 中 ,qd,(i-s,i+s) 为 降 维 后 图 像 4;., 与 w; 的 欧 
氏 距 离 , 即 qd,(i- s,i - s) = || os. -4,6 为 原始 图 像 u 在 低 
维 流 形 中 的 映射 ,w, 为 权重 值 , 取 w, =e ”22 oy 
7 为 正 整 数 , 取 m = round ie 

归 一 化 公式 (4) ,得 到 最 终 的 距离 函数 表达 式 : 

max(d)- min(d) 

在 一 个 心动 周期 中 , 当 w 处 于 心脏 舒张 未 期 时 , 距 
离 函 数值 达到 最 小 。 图 4 为 一 组 长 度 为 1024 的 IVUS 
图 像 序列 距离 函数 曲线 ,在 一 定 程 度 上 反映 了 心脏 周期 
性 的 运动 规律 。 
1.2.4 提取 关键 帧 距离 函数 的 局 部 最 小 值 对 应 序列 图 
像 关 键 帆 的 位 置 ,本 文采 用 以 下 方法 寻找 距离 函数 局 部 
最 小 值 : 

(1) 搜 索 起 始点 。 以 距离 函数 4 的 最 小 值 作为 搜索 


型 ,频率 40 MHz, 回 撤 速 度 0.5 
mm/s, 采 样 频率 为 30 帆 /s。 实 验 
SN 平台 采用 个 人 计算 机 (CPU: 英 
FEAR SB — XS TE 13-3240 @ 
3.40 GHz 双核 ,4 GAF), 


gaga 

0.02. 0.040.060.08 Matlab2011a 编 程 。 

4-0.020 9. D 2.1 IVUS 图 像 序列 参数 定量 
分 析 与 比较 


临床 诊断 和 治疗 冠 心 病 ， 
第 用 到 的 诊断 指标 中 ,血管 容 
TR EPIS CER, XE fur TEUER SU LAT HRS 
数 ,这 些 参数 的 定量 计算 需要 检测 外 弹力 膜 (external 
elastic membrane, EEM) 界面 和 管 腔 内 膜 
(lumen-intima) 界 面 ,EEM 界面 近似 为 血管 壁 , 管 腔 内 
膜 界面 近似 为 管 腔 边 界 。 但 目前 未 存在 精准 的 提取 这 
两 个 界面 的 方法 ,在 临床 医生 的 指导 下 ,人 工 标 记 这 两 
个 界面 。EEM 界 面 可 近似 为 椭圆 ,在 IVUS 图 像 中 标 
记 EEM 的 八 个 边界 点 ,最 小 二 乘 拟 合 出 椭圆 的 一 般 方 
程 ,以 确定 血管 面积 及 血管 最 大 最 小 直径 。 管 腔 内 膜 界 
面 不 能 用 确定 的 模型 模拟 ,人 工 标 记 管 腔 边界 。 管 腔 中 
DELH”: 
0, -[i, v] 7 x unva (6) 
公式 (6) 中 ,[j,vi 为 边界 点 的 坐标 ,NN 为 管 腔 / 内 膜 
边界 点 数 , 以 此 来 确定 管 腔 面积 .最 大 最 小 直径 。 结 合 
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ca 序列 的 血管 面积 , 管 腔 面积 变化 曲线 。 红 色 实 线 为 原始 
0.8 


0 200 400 600 800 1000 
frames 


图 4 一 组 IVUS 图 像 序列 距离 函数 曲线 


Fig.4 Distance function of an IVUS sequence. 


两 界面 参数 获得 斑 块 负荷 等 血管 参数 。 原 始 序列 和 门 
控 序列 血管 容积 分 别 定义 为 : 

Vov 20.5d X 17 (Soy) + Sov(i+t))) (7) 

Vav 20.54 X Y hv À+ Sevi + D) (8) 

公式 (7,8) 中 ,LL,! 分 别 是 原始 序列 和 门 控 序列 的 帧 
ZE, Sov 0) ,SovQ 人 分 别 是 原始 序列 和 门 控 序列 第 帧 图 像 
的 血管 横 截 面积,@ 为 一 帧 图 像 对 应 的 血管 厚度 ,本文 
d= 1/60(mm),hQ) 为 相 邻 门 控 帧 之 间 的 帧 距 。 原 始 序 列 
和 门 探 序列 管 腔 容 积 的 计算 类 似 。 

斑 块 负 答 定义 为 : 

B,=(S,-S)/S, (9) 

公式 (9) 中 ,5, 为 血管 横 截 面积 ,S 为 管 腔 横 截面 积 。 

在 临床 医生 的 指导 下 ,人 工 勾 画 13 组 IVUS 图 像 序 
列 的 EEM 界 面 和 管 腔 内 膜 界面 。 计 算 门 控 前 后 序列 的 
血管 容积 , 管 腔 容 积 . 平 均 斑 块 负荷 ;检验 原始 序列 和 门 
控 序 列 这 三 个 参数 的 差异 显著 性 ;对 每 一 组 IVUS 原始 
序列 和 门 控 序列 的 血管 横 截 面积 方差 , 管 腔 横 截面 积 方 
差 进行 右边 检验 。 参 数 差异 显著 性 检验 采用 配对 秩 和 
检验 方法 。 原 始 和 门 控 序 列 血 管 横 截 面积 方差 右边 检 
验方 法 是 : 设 原始 序列 和 门 控 序列 的 血管 横 截 面积 序列 
分 别 为 X,Y ,方差 0,62。 做 出 检验 假设 , 原 假设 Ho: 
01> 052, 备 择 假设 Hi:01 友 63。 先 对 和 ,了 的 分 布 进行 k 检 
验 ,检验 它们 是 否 符合 正 态 分 布 ; 若 X ,7 都 符合 正 态 分 
布 ,再 对 这 两 个 正 态 分 布 的 方差 进行 右边 检验 (F 检 
验 )。 管 腔 横 截 面积 方差 的 右边 检验 类 似 。 实 验 统计 结 
果 如 表 1 所 示 , 门 控 序 列 的 血管 容积 , 管 腔 容 积 显 著 小 
于 原始 序列 。 同 时 , 门 控 前 后 的 平均 斑 块 负荷 无 显著 性 
差异 ,说 明 门 控 序 列 保留 了 原始 序列 中 管 腔 面 积 和 血管 
面积 的 相关 信息 ,符合 临床 诊断 需要 。 并 且 , 通 过 原始 
序列 和 门 控 序 列 的 血管 横 截面 积 方差 管 腔 横 截 面积 的 
方差 进行 右边 检验 得 出 : 门 控 序 列 的 血管 面积 方差 和 管 
腔 面积 方差 小 于 原始 序列 , 即 门 控 序列 的 管 腔 面 积 及 血管 
面积 变化 比 原始 序列 小 ,表明 门 控 序 列 比 原始 序列 稳定 。 


序列 血管 面积 变化 曲线 ,红色 虚线 为 原始 序列 管 腔 面积 
变化 曲线 , 蓝 色 实 线 为 门 控 序列 血管 面积 变化 曲线 , 蓝 
色 虚 线 为 门 控 序列 管 腔 面积 变化 曲线 , 蓝 色 星 号 对 应 关 
键 帧 血管 , 管 腔 面积 。 门 探 序列 面积 变 比 原始 序列 平缓 ， 
保留 了 原始 序列 中 的 信息 ,具有 更 加 直观 的 视觉 效果 。 

综 上 所 述 ,本 文 方法 提取 到 的 关键 帧 组 成 的 门 控 序 
列 的 血管 容积 管 腔 容积 显著 小 于 原始 序列 , 门 控 前 后 
序列 的 平均 斑 块 负荷 差异 性 不 显著 ,满足 临床 诊断 要 
求 。 门 控 序 列 的 血管 面积 方差 和 管 腔 面积 的 方差 显 车 
小 于 原始 序列 ,表明 门 控 序列 较 原始 序列 稳定 。 
22 IVUS 图 像 序列 纵 切 图 

图 6 为 一 组 IVUS 序列 在 不 同 的 纵 切 方向 上 (0。， 
45°,90°,145°) 的 纵 切 图 ,上 部 为 原始 序列 ,下 部 为 门 控 
序列 。 图 7 为 四 组 序列 关键 帧 检索 的 结果 ,上 部 为 原始 
序列 ,下 部 为 门 控 序列 。 从 纵 切 图 上 观察 到 ;由 关键 帧 
所 组 成 的 门 控 序 列 的 管 壁 . 管 腔 、 斑 块 与 原始 序列 具有 
高 度 的 相似 性 ,血管 内 腔 和 管 壁 轮廓 和 斑 块 比 原始 序列 
要 平滑 ,连续 性 更 好 ,抑制 了 锯齿 状 振 荡 ,具有 更 好 的 视 
2.3 与 已 有 方法 比较 

文献 [6] 提 出 一 种 基于 IVUS 图 像 序列 灰 度 的 门 控 
方法 ,其 3 个 步骤 为 :(1) 选 择 每 帧 图 像 像素 点 局 部 均值 
为 参考 点 , 沿 序列 形成 一 维 信号 复 ,通过 幅 值 滤波 和 光 
流 滤波 保留 与 心脏 运动 信息 相关 性 很 强 的 一 维 信号 簇 ; 
(2) 带 通 滤 波 一 维 信号 艇 ,组 成 一 个 心动 周期 剖面 图 ; 
(3) 寻 找 滤波 后 的 一 维 信和 号 徐 极 值 ,提取 关 键 帧 。 该 方 
法 能 快速 准确 的 提取 门 控 序列 。 本 文 方法 与 文献 [6] 所 
提出 的 方法 进行 比较 ,和 N 为 实验 序列 长 度 ,分 别 统计 了 
两 种 方法 得 到 的 门 控 序列 中 关键 帧 帧 数 (N,、Nia) ,关键 
帧 平均 距离 Ey, 以 及 关键 帧 之 间 的 距离 占 心脏 周期 
的 比例 分 数 EF, ,结果 见 表 2, 可见 两 种 方法 得 到 的 门 
控 序 列 中 的 总 帧 数 接近 ,Ey,EF, 值 较 小 ,说 明 两 种 方 
法 检索 到 的 关键 帧 比较 一 致 。 但 文献 [6] 中 方法 对 华 通 
滤波 器 参数 敏感 ,难以 提取 心率 不 齐 患者 及 冠 脉 血管 结 
构 多 变 的 IVUS 图 像 序列 的 门 控 帧 。 


3 结论 

为 门 控 序 列 包 含 了 原始 序列 中 的 关键 信息 ,提取 
关键 帧 门 控 序列 ,降低 了 临床 血管 数据 的 分 析 量 , 相 比 
于 原始 序列 , 门 探 序列 具有 更 好 的 视觉 效果 。 本 文 提出 
了 一 种 简单 的 ,稳定 的 基于 IVUS 图 像 序列 的 门 控 方 
法 。 应 用 流 形 学 习 的 方法 ,将 高 维 的 IVUS 数据 降 到 低 
维 流 形 中 ,减少 了 计算 机 存储 和 分 析 量 。 以 低 维 流 形 中 
的 特征 分 量 作 为 显著 性 特征 ,构造 距离 函数 ,检索 并 提 
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表 1 13 组 IVUS 图像 原始 序列 及 关键 帧 门 控 序列 的 血管 参数 和 分 析 统 计 
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Tab.1 Vascular parameters and statistic analysis on 13 IVUS non-gated versus key-frames gated sequences 


No Segment length Vessel volume(mm?) Lumen volume(mm ) Mean plaque burden(%) Variance test of Variance test 

(frames) Original Gated Original Gated Original Gated vessel area of lumen area 
No.1 891 175.87 174.41 146.56 143.95 16.66 17.46 HO HO 
No.2 883 268.37 265.63 239.08 236.50 10.91 10.97 HO HO 
No.3 916 131.89 123.62 118.41 109.80 10.22 11.18 HO HO 
No.4 711 194.96 191.71 94.83 94.82 51.36 50.54 HO HO 
No.5 912 107.28 103.93 50.61 48.85 52.82 53.00 HO HO 
No.6 873 430.72 427.64 270.01 264.92 37.31 38.05 HO HO 
No.7 915 133.97 129.66 94.10 90.67 29.76 30.07 HO HO 
No.8 893 121.52 119.88 103.69 102.18 14.68 14.77 HO HO 
No.9 785 120.41 119.45 81.82 81.25 32.05 31.98 HO HO 
No.10 933 196.49 191.98 117.21 117.86 40.35 38.61 HO HO 
No.11 890 303.50 296.11 273.277 267.41 9.96 9.69 HO HO 
No.12 921 256.41 246.32 215.28 206.88 16.04 16.01 HO HO 
No.13 748 107.29 103.83 52.89 51.94 50.70 49.98 HO HO 
Significant difference test 0.001 0.003 NS 


The confidence level of all hypothesis test is 0.05. NS indicates no significant difference between the original and the gated sequences. HO represents 


a greater variance of the vessel cross-sectional area or the variance of the lumen cross-sectional area of the original sequence than the gated ones, 


while H1 indicates the opposite conclusion. 
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图 5 原始 序列 和 关键 帧 门 控 序 列 血管 及 管 腔 面 积 变化 曲线 


Fig.5 Vessel and lumen cross-sectional area of an original 


and key-frames gated set. 


取 心 脏 舒张 末期 的 IVUS 图 像 ,组 成 关键 帧 门 控 序列 。 
本 文 算法 的 优势 体现 为 两 方面 :(1) 存 储 和 计算 量 减少 ; 


(2) 对 于 心率 不 齐 的 患 
形 学 习 中 的 特征 分 量 含 有 复杂 的 运动 规 和 


,能 准确 地 寻找 门 控 序 列 。 流 


,如 图 像 局 部 


平移 ,旋转 .尺寸 变化 等 运动 参数 ,本 文 仅 针 对 IVUS 图 
像 序列 纵向 运动 伪 影 的 抑制 进行 研究 ,对 短 轴 旋转 伪 影 


的 抑制 和 补偿 有 竺 进一步 探讨 。 
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Fig.6 An original IVUS set longitudinal cuts at different angles (0°, 45°, 90°, 145°) and its gated cuts. 


图 7 4 组 IVUS 原始 序列 及 其 门 控 序 列 的 纵 切 图 
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Fig.7 Four original IVUS sequences and their longitudinal cuts and gated cuts. 


R2 13 本 文 方法 与 文献 [6] 方 法 对 比 结果 


proposed in reference" 


Tab.2 Comparative results of our method and the method 
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